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Abstract. Mode identification is crucial for an asteroseismological study of 
any significance. Contrarily to spectroscopic techniques, methods such as period- 
fitting and multi-colour photometry do not provide a full reconstruction of non- 
radial pulsations. We present a new method of spectroscopic mode identification 
and test it on time-series of synthetic spectra appropriate for pulsating subdwarf- 
B stars. We then apply it to the newly discovered slowly pulsating subdwarf-B 
star HD 4539. 



1. Introduction 

Subdwarf B (sdB) stars are low-mass (~ CSA/q) core helium burning objects 
with thin and mostly inert hydrogen-ri ch residual e nvelopes. They belong to 
the Extreme Horizontal Branch (EHB) (lHebeij[l986l ) and remain h ot (20000 < 



Teff < 40000) and compact (5 < log g < 7) throughout their lifetime (jSaffer et al 



19941 ). before evolving towards the white dwarf cooling sequence without expe- 
riencing t he Asymptotic Giant B ranch and Planetary Nebula phases of stellar 
evolution. iHan et al.1 (|2002l . l2003l ^ use binary population synthesis calculations 
to demonstrate the formation of sdB stars through several possib le channels, but 
result ing models require furth er comparison with observation. [Kilkenny et al 



( 19971 ) and lGreen et al.l (j2003l ) respectively discovered that some sdB stars un- 
dergo fast and slow non-radial pulsations. This means asteroseismology can be 
used to study the internal structure of these stars and so constrain evolution 
models. 

After having measuring the pulsation periods, the first step of any seismic 
modelling is to identify the pulsation modes that are excitecfl, which is far from 
being trivial. In the case of pulsating sdB stars, most m ode identifications up 
to date have been pe r formed using either p eriod fitting ([Brassard et al.1 l200ll : 
Charpinet eraP l2005l: iRandall et all l2006allbl. etc.) or multicolour amplitude 
ratios (|Jefferv et al.ll2004l . l2006l . etc.). However, these methods do not provide a 
full reconstruction of the pulsations, and should therefore be complemented by 
spectroscopic mode identification techniques, such as the one we present here. 



^ Since non-radiai puisations are modelied by splierical harmonics l^^m, mode identification 
consists in assigning values to the sphericai wavenumbers n, i and m, respectiveiy the number 
of nodes of the radial dispiacement, the number of nodai fines on the steffar surface and the 
number of such fines passing through the rotation axis of the star. 
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2. Spectroscopic Mode Identification 



As reported earlier ( Schoenaers &: Lvnas-Gravll2006l ) . line-profile variations (Ipv), 



although weak, are to be expected in pulsating subdwarf-B stars. Furthermore, 
they have been observed, for instance in PG 1325+101 bv iTelting fc 0stensenl 
([2004). A reliable spectroscopic mode identification method is therefore needed 
to take advantage of these and similar observations. 

2.1. Why not use the moment method? 

Nowadays, the most commonly used spectroscopic mode identification method 
i s the "moment method" (hereafter MM) . It was first introduced by iBalona 



( 1986al lbl. Il987l ) and later reformulat ed into computer-f riendl y terms and ap- 



phed to Main Sequence p ulsators bv lAerts et al.l ( 19921 ) and lAertsI (|l996l ) for 



monoperiodic stars, while iBriquet &: AertJ l 20031 ) extended the formalism to 



multiperiodic stars. In this method, one does not try to identify pulsation modes 
directly from the line-profiles, because this is often impossible for multiperiodic 
stars, but instead replaces each observed line-profile by its first few moments 
and compares their time dependence with that of theoretical moments. The 
latter depending on the wavenumbers i and m of the pulsation mode, as well 
as on the pulsation amplitude and on the inclination angle i of the star, this 
method should in principle allow a complete reconstruction of any non-radial 
pulsation. It has been us ed successfully for the study of Cephei and Slowly 
Pulsating B (SPB) stars (|Briauet AertsI l2003l : fPe Cat et al.ll2005l 'l. However, 



this method makes use of a crucial (and not always justified) approximation: 
except for radial velocity variations, the MM neglects all other variations due to 
non-radial pulsations. 

However, because of non-radial pulsations, different points on the stellar sur- 
face not only have different radial velocities, but also different temperatures, log 5 
and orientations, and their contributions to a line-profile therefore have different 
amplitudes. Hence, Ipv do not only come from the Doppler shift in wavelength, 
but also from the complex temperature and log g behavior on the stellar surface. 
This is especially the case in sdB pulsators: because they are so dense, any signif- 
icant displacement of the stellar surface due to pulsation is kept at its minimum, 
and instead the available energy causes temperature perturbations in the stellar 
photosphere. To demonstrate that neglecting temperature, gravity and geomet- 
ric variations can lead to a mistaken mode identification, even in the simplest 
cases, we used the MM to identify various monoper iodic pulsations present in 
time-s eries of synthetic spectra computed following Schoenaers Sz Lvnas-Gravl 



(|2006l 1. On the' one hand, it could be seen that the MM always provide, within 
its fiv^ best solutions, the proper mode identification (as well as very good es- 
timates of i, Ap and Veq) when only radial velocity variations were taken into 
account. However, on the other hand, when temperature, gravity and geometric 
variations were taken into account (as indeed they should be) this test demon- 
strated that even in the simple case of monoperiodic pulsations, the current MM 
falls short of properly identifying all pulsational characteristics. Clearly, in the 
case of pulsating sdB stars, where temperature, pressure and geometric effects 



^Following iBriguet fc AertsI (|2003f l's advice. 
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have to be taken into account, an improved version of the moment method is 
needed. 

2.2. The "synthetic moment" method 

Schoenaers &: Lvnas-Gravl ( 20061 ) introduced a state-of-the-art code that com- 



putes times-series of synthetic spectra and Ipv. The "synthetic moment" method 
(hereafter SMM) introduced here takes advantage of this abihty to compute 
"synthetic moments" for these extremely accurate synthetic Ipv, and to com- 
pare them to observed moments: for a given target and observed hne-profile 
variation, a grid of time-series of synthetic spectra spanning a suitable parame- 
ter space is computed, the first three moments of each of these synthetic Ipv are 
obtained, and the best fit between observed and synthetic moments is found by 
minimising a merit function of the form 

+ I < >syn,j - < >obs,j I (1) 
+ ( < f ^ >syn,j - < f ^ >obs,j 




where A'obs is the number of observations, odd moments have average zero, and 

*syn,i 

,,2 



the synthetic second moment < v'^ >syn,i is adjusted to have the same average 



as its observed counterpart < v >obs,i- 

The choice of the parameter space to be searched is crucial: if too limited 
one might miss the proper mode identification, but if too broad the computing 
time of the synthetic Ipv can become prohibitive. One could consider building 
an extensive grid of synthetic Ipv by varying these parameters more or less "con- 
tinuously", but because most parameters are target-dependent, it is simply more 
efficient to compute one grid for each pulsating star. Provided spectroscopic ob- 
servations of the target are obtained at high enough SNR, Teg, log 5, Ueqsinz and 
[X/H] can be determined to good accuracy. Before attempting mode identifica- 
tion, one should also obtain reliable values of the pulsation period(s) P, either 
through (ideally simultaneous) photometry or radial velocity measurements. If 
the latter are obtained (either directly through spectroscopic measurements or 
indirectly from the first moment), then good estimates of the physical velocity 
amplitudes of the pulsation mode(s) can be obtained; if not, then pulsation am- 
plitudes can in principle remain a free parameter, along with the inclination i' 
and the actual mode identification {i',m'). 

Ideally, when dealing with a multiperiodic star, all pulsation modes should 
be identified simultaneously, but this would require a very large grid of synthetic 
Ipv, which can become prohibitive. It is however possible to attempt mode 
identification as follows: 

1. Use the SMM on the dominant pulsation mode with a suitable primary grid 
of synthetic Ipv to obtain , m'l ) and i' for the five best fitting solutions 
(minimising A^/ „j/ ). Then, if at all possible, further discriminate between 
these five best solutions (see next subsection for more details). 
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2. For each of the possible i' and {£[,m[) combinations obtained in Step 1, 
build a secondary grid of synthetic Ipv (one for each {£'2, 1712) combination) 
and use the SMM to identify the second pulsation mode. 



3. Step 2 can be repeated at will until all modes have been succesfully iden- 
tified. 



2.3. Discriminating between the few best solutions 

As experience ha s show n in the case of Main Sequence pulsators (see for exam- 
ple Ee Cat et al.l it is in general safer not to accept straight away the 
best solution of the MM (and by extension, of the SMM) as the proper mode 
identification, but instead to discriminate between a few best solutions using 
some other mode identification method. Similarly, to confirm the identification 
of pulsation modes in PG1716 stars, we suggest a visual inspection of the fit 
between observed and synthetic moment s, as well as the use of the IPS (for 
Intens ity Period Search) diag:nostic (Schriivers et al.l Il997l : iTelting &: Schriivers 



ity^ . ^ 

1997al lbl: ISchriivers k Telth^ Il999l ^ to compare the amplitude and phase dia- 



grams of observed Ipv with those of synthetic profiles corresponding to the few 
best solutions of the SMM. 



2.4. Test on synthetic data 

Figure [U shows in grey-scale a biperiodic Ipv, representative of the PG1716 star 
PG 1338+481 (Teff = 27 000 K, logg = 5.4, M = 0.48 Mq, R = 0.23 Rq, Veq = 
5km • s~^, [He/H]= —1.3, [Mg/H]= —1), that was used (amongst many others) 
to test the SMM. Pulsation characteristics are given in the caption. Appropriate 





Figure 1. Grey-scale plots of synthetic Ipv {left: He I A4471 A and Mg II 
A448lA, right: He I A6678A) with periods and i'-values representative of 
PG 1338+481 l|Randall et al.ll2006al) : i = 60°, (4,toi) = (1,1), Pi = 3828s, 
Ap,i = 2.0km-s-i and (^2,^2) = (1,0), P2 = 3530s, Ap,2 = 1.4km- s^^. 
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non-adiabatic parameters R and tpx were computed using iRandall et aD (|2005l Vs 
Equations 46 and 470 

Applying the SMM to these synthetic Ipv of He I A6678 ^ provides, for the 
first pulsation mode, the five best solutions listed in the first four columns of 
Table [TJ Visual inspection of the match between "observed" and synthetic mo- 



Table 1. Five best identifications obtained by the SMM for the biperiodic 
pulsation modelled above. The first column is the inclination of the rotation 
axis with regard to the line-of-sight (in degrees), the second and third columns 
give the {£[,m[)-va\ues of the first mode, and the fourth provides the value 
of Af'^ m'^, minimum for the best fitting solution. Columns 5 to 7 give the 
{£'2, m2)-values and the value of A^^ for the second mode, provided the first 
was properly identified. For each mode, the proper identification is written in 
bold font. 



First mode Second mode 



i' 










1^2 




60 


1 


1 


4.5436200 


1 





1.31584586 


30 


2 


1 


5.1651505 


1 


-1 


1.90534104 


30 


3 


1 


5.2640002 


3 


2 


2.25563231 


85 


2 





5.4004608 


3 


-3 


2.43201440 


60 


3 





5.4952073 


3 


-2 


2.44369818 



ments (not plotted here for the sake of brevity, but see lSchoenaer i ([2008l Vs Figure 
4.23) allowed us to rule out all but the {i'i,mi) = (1, 1) and (3, 1) identifications 
for the first mode. Further discrimination can be achieved by comparing the 
"observed" IPS diagnostic (shown in Figure [2D to those of the {i'i,m'i) = (1, 1) 
and (3,1) identifications (plotted in Figure [3]). The phase diagrams of (3,1) do 
not exhibit the right behaviour, while those of (1,1) do. The latter is therefore 
selected as the identification of the dominant mode. 

Keeping (£i,mi) = (1,1) and i = 60° fixed, and applying the SMM to 
identify the second pulsation mode yields the five best identifications listed 
in Columns 5 to 7 of Table [H The IPS diagnostic of the best identification 
{£'2,1712) = (1,0) plotted in Figure H] matches the "observed" IPS diagnostic very 
well, and this shows the SMM does provide a reliable mode identification, even 
for multiperiodic stars. 



^Note tliat Equation 47 of lRandalt et at.l (120051 ) sliouid read (iPt) = 7r-2.30exp(-P^'''^**/1800). 

^ Simitar resufts ar e obtained wfien appiying tlie SMM to He I A4471 A and Mg II A4481 A 
(|Schoenaers|[200l ). 
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Figure 2. "Observed" IPS diagnostic for the Ipv shown in Figure [T] The 
two leftmost panels show the amplitude (left) and phase (right) of variations 
with the main frequency across the profile, while in the two rightmost panels 
the amplitude (left) and phase (right) plots are for variations with the first 
harmonics of the main frequency. 
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Figure 3. Same as Figure [H but for possible identifications of the first 
pulsation mode: (i[,m[) = (1,1) (top) and (3,1) (bottom). 
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Figure 4. Same as Figure [21 but for possible identifications of the first 
pulsation mode: (£[,m[) — (1,1) (top) and (3,1) (bottom). 
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Figure 5. Line- profile variations in HD 4539 observed in August 
2005 (two leftmost panels) on MJD=53611. 13918 {solid line) and 
MJD=53612. 03206 {dashed line), and in August 2006 {two rightmost panels) 
on MJD=53958. 15420 {solid line) and MJD==53952. 16671 {dashed line). 



3. Application to a New Slowly Pulsating Subdwarf-B Star: HD 4539 

HD 4539 {B ~ 10.12, 2000.0 coordinates: a = 0^^47^29.2^, 5 = + 09° 58^56'^) is 



one of the brightest kno wn sdB stars. Its surface parameters (see e.g.lBascheck et al 



97li 



Saffer et al.lll99J ) are similar to those of PG1716 stars, and Schoenaers k, Lynas-Gray 



( 20071 ) report the discovery of non-radial pulsations in HD 4539, together with 
their frequency analysis, but without mode identification and subsequent anal- 
ysis. In Figure [5l we plot Ipv observed in HD 4539 in August 2005 and 2006. 
These Ipv are much stronger than the average noise present in our spectra, indi- 
cating they are r eal. A careful radial velocity ana lysis and period search yielded 



four frequencies ( Schoenaers &: Lvnas-Gravll2007l ). of which we identify the dom- 



inant one (Pi = 9310 s) in this preliminary report. A more thorough analysis 
will be presented elsewhere. 

The wavelength range of our spectra includes He I A 4471 A and Mg II 
A 4481 A, and we therefore applied the SMM to the first three moments of 
these two lines. In Table [2] are listed the five best-fitting identifications of the 
dominant pulsation mode. These five best solutions (including the inclination) 



Table 2. Five best fitting SMM identifications of Pi in HD 4539, using He 1 
A 4471 A {left) and Mg II A 4481 A {right). In each case, the first and second 
columns give the tentative {£[,m[) identification, the third column provides 
the corresponding inclination, and the fourth is the merit function A^'^ „i'^. 





He I A 4471 A 




Mg 
m'^ 


II A 4481 A 


1 


1 


85 


35.7424808 


3 


3 


85 


55.5834533 


2 


1 


60 


35.7468680 


3 


1 


30 


55.5841337 


2 





85 


35.7529125 


2 





85 


55.5845978 


3 


3 


85 


35.8536815 


2 


1 


60 


55.5993271 


3 


1 


30 


35.8943103 


1 


1 


85 


55.6007015 
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are the same (although found in a different order) whether the identification is 
performed using He I A 4471 A or Mg II A 4481 A, which aheady indicates they 
are quite robust. Visual inspection of the moments did not allow us to discard 
any particular identification, and therefore we compared the IPS diagnostic of all 
five best-fitting solutions to the observed IPS diagnostic of HD 4539, computed 
for both He I A 4471 A and Mg II A 4481 A. The same conclusion could be 
drawn in both cases, namely that the dominant pulsation mode (Pi = 9310 s) 
in HD 4539 can be identified as a = 1, rrii = 1 mode, while the inclination of 
the rotation axis with respect to the line-of-sight is (close to) 85°. 
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